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Planck-2G15 data seem to favour a large value of the lensing amplitude parameter, A-l = 
1.22 ± 0.10, in CMB spectra. This result is in 2a tension with the lensing reconstruction result, 
A’^ = 0.95 ± 0.04. In this paper, we simulate several CMB anisotropy and CMB lensing spectra 
based on Planck-2015 best-fit cosmological parameter values and Planck blue book beam and noise 
specifications. We analyse several modified gravity models within the effective field theory frame¬ 
work against these simulations and find that models whose effective Newton constant is enhanced 
can modulate the CMB anisotropy spectra in a way similar to that of the j4l parameter. However, 
in order to lens the CMB anisotropies sufficiently, like in the P/ancfc-2015 results, the growth of 
matter perturbations is substantially enhanced and gives a high as value. This in turn proves to be 
problematic when combining these data to other probes, like weak lensing from CFHTLenS, that 
favour a smaller amplitude of matter fluctuations. 


I. INTRODUCTION 

Based on the full-mission Planck observations of tem¬ 
perature and polarization anisotropies of the cosmic mi¬ 
crowave background (CMB) radiation, Planck-2015 re¬ 
sults show that the temperature and polarization power 
spectra are consistent with the standard spatially-flat six- 
parameter ACDM cosmology with a primordial power- 
law spectrum of adiabatic scalar perturbations. Here¬ 
after we shall call this model the base-ACDM. On the 
other hand, the same data, especially the temperature- 
temperature (TT) spectrum reveals some tension with 
the CMB lensing deflection angle (d) spectrum recon¬ 
structed from the same maps. In details, the lensing 
amplitude in CMB temperature and polarization spec¬ 
tra, Al = 1.22 ± 0.10, is in 2a tension with the am¬ 
plitude of the CMB trispectrum reconstructed lensing 
deflection angle spectrum, Aff = 0.95 ± 0.04 while it 
is expected that in the base-ACDM model both these 
quantities should be equal to unity. 

The Planck collaboration finds that, compared with 
the base-ACDM model, the base-ACDM-|-AL model can 
reduce the logarithmic likelihood (Ay^ = —6.1) and pro¬ 
vide a better fit to the data sets with ^ Al = 1.28 or 
marginalized constraint Al = 1.22 ± 0.10 [T] . More im¬ 
portantly, they find that there is roughly equal preference 
for high Al from intermediate and high multipoles (i.e., 
the Plik likelihood; Ay^ = —2.6) and from the low-f like¬ 
lihood (Ay^ = —3.1) with a further small change coming 
from the priors. This means that the base-ACDM-pAL 
model can provide a better fit than base-ACDM model 
against both TT and lowP data sets. However, the in- 


^ From the PLA-PR2-2015 official chains 

base-Alens-plikHM-TT-lowTEB at http://pla.esac.esa. 
int/pla/ 


crease in Al will induce changes on the full sets of cosmo¬ 
logical parameters as mentioned in the reference [T]. For 
example, compared with the base-ACDM fit, the scalar 
index, Ug, is increased by 1%, the primordial scalar spec¬ 
trum amplitude, Ag, is reduced by 4% and the effective 
amplitude of the TT spectrum, AgC”^"^, is reduced by 
1%. Through the complicated relationship between pa¬ 
rameters and their degeneracy the re-ionization optical 
depth parameter, t, falls to 0.060, which is roughly in 2cr 
tension with Planck-2013 temperature -I- WMAP low-.^ 
polarization data results of r = 0.089lg;Q(4. 


Inspired by these observations, in this paper, we inves¬ 
tigate whether some modifications of gravity can relieve 
the tension between Planck CMB anisotropy spectra and 
CMB lensing results. To do so we will simulate a tension 
CMB data set that resembles the tension present in the 
Planck-2015 results and we will try to fit the resulting 
power spectra with different models to have a glimpse of 
the changes in the parameter that arise because of this 
tension. 


We shall show that the modified gravity models that we 
consider are not capable of alleviating this kind of tension 
without affecting substantially the growth of structure. 
These models will then have problems in fitting simul¬ 
taneously other cosmological probes. Moreover we will 
show that, if the amplitude of the reconstructed CMB 
lensing power spectrum is in tension with the amount 
of CMB lensing in the anisotropy spectra, these kind of 
models struggle to relieve it. As we will elaborate later 
this happens because modified gravity models predict the 
same amplitude for the reconstructed CMB lensing po¬ 
tential and the lensing effect on the CMB so that the ten¬ 
sion is usually moved to other parameters of the model. 
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II. MOCK DATA AND FIDUCIAL 
PARAMETERS 


We analyse several modified gravity models against 
two sets of simulations of CMB spectra (TT,TE,EE) and 
CMB lensing spectra (dd and Td) with the fiducial cos¬ 
mological parameter equal to the Planck-2Qlh data re¬ 
lease best-fit values and Planck blue book beam and noise 
specifications. The motivation of doing these simulations 
are mainly two. Eirst of all the Planck-201b likelihood 
code and the corresponding spectrum data are not yet 
publicly available. Second, because the cosmological pa¬ 
rameters are degenerate with each other in a complicated 
way, by using simulations, we can efficiently isolate and 
study the effects coming from different parameters and 
their combinations. 

In the following section, we will briefly review the Eu- 
turCMB^ package [2] that we used to produce our simu¬ 
lations. The construction of the power spectrum simula¬ 
tions are mainly made by two essential parts, firstly the 
mock data generation and the likelihood construction. 

Eor the mock data generation we assume, for simplic¬ 
ity, the raw maps to be composed of Gaussian signal 
(s^m) and a uniform Gaussian white noise (n^m)) 

= ^£m + ^fm i (1) 

where the superscript P stands for temperature and E- 
mode polarization as we do not consider B-mode polar¬ 
ization. The observation phase then convolves this sig¬ 
nal with a Gaussian beam characterized by the full width 
half maximum parameter dfwhm- On the top of it, we also 
need to add a Gaussian white noise. Up to a normaliza¬ 
tion, the noise power spectrum can be approximated as 




1 ) 


^fwhm 

8 log 2 


( 2 ) 

where ap models the root mean square of the instrumen¬ 
tal noise. In this paper, we shall adopt the Planck blue- 
book [3] beam and noise parameters. The exact values 
used in our simulations are listed in Table [D The mock 
data generator package EuturCMB automatically com¬ 
putes the noise power spectra of the lensing deflection 
angle based on the Hu-Okamoto [1] quadratic estimator 
algorithms from the given temperature and polarization 
spectra and noises. 

Then, we feed the fiducial cosmological model into Eu¬ 
turCMB. In this paper, we generate the fiducial spec¬ 
tra Cj^, Cf^, as well as Cj'^ from the pub¬ 

lic Boltzmann code CAMB 0 [g. In order to mimic 
Planck-2015 results, the cosmological parameter used to 
produce the fiducial spectra are set to the best-fit val¬ 
ues of the base-ACDM-|-AL model to the TT-|-lowP and 


^ http://lpsc.in2p3.fr/p6rotto/ 


TABLE I: Planck blue book instrumental specifications 


Experiment 

Frequency 

^beam 

(Tt 

(Jp 


Planck. 

217 

5.02 

13.1 

26.7 



143 

7.30 

6.0 

11.4 



100 

9.68 

6.8 

10.9 



Frequencies in GHz. Beam size 0beam is the FWHM in arcmin- 
utes. Sensitivities ap and ap are in fiK per FWHM beam. 


TT-|-lowP-|-lensing data set of the Planck-2015 data re¬ 
lease (see Table|^. With these input the EuturCMB code 
will generate Gaussian-distributed random fields for 
and estimate the corresponding spectra ■ For fur¬ 
ther details about the numerical implementation of the 
EuturCMB code we refer the reader to [5]. 

The second step is to construct the likelihood for the 
mock data sets. Since it is generated from a Gaussian 
realization, we can write the likelihood function of the 
data given the theoretical template as [7] 

/:(a|0) oc exp [C(©)^^]a^ , (3) 

where a = *^ata vector, 0 = 

(01,02, ■ • •) is the parameter vector, and (7(0) is the the¬ 
oretical data covariance matrix. For the detailed con¬ 
struction of the spectrum likelihood we refer to HIS]. 

After the above operations, we build two mock data 
sets (Cj^, Cf^ and that we shall 

call Mock-A and Mock-B, whose fiducial cosmological pa¬ 
rameter values (see Table ITT]) are, respectively, the best- 
fit values of base-ACDM-I^L to Planck-2015 TT-|-lowP 
and Planck-2015 TT-blowP-blensing data sets^. Since the 
Mock-A data set, which mimics Planck-2015 TT-|-lowP, 
is generated from Ap ~ 1.3, we can treat it as a realiza¬ 
tion of a non-AGDM universe; Mock-B data, that mimics 
Planck-2015 TT-|-lowP-flensing, Ap ^ 1.0, is closer to a 
realization of ACDM universe. 

Based on Mock-A and Mock-B data sets, we build a 
“tension” data set, called Mock-C, by combining ((7j^, 
Cf^, Cf^) from Mock-A and {Cf‘^ and Cj'^) from Mock- 
B. The resulting data set should mimic the data compila¬ 
tion of TT-|-lowP-|-lensing data sets of Planck-2015 data 
release. In the rest of the paper, we will study several 
modified gravity models against Mock-A and Mock-C 
data sets to see whether the modified gravity models can 
or cannot reconcile the tension between CMB anisotropy 
data and CMB lensing data. 


® The best-fit values are read from the PLA-PR2- 
2015 official chains base-Alens-plikHM-TT-lowTEB and 
base-Alens-plikHM-TT-lowTEB-lensing at http://pla.esac. 
esa.int/pla/ 
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FIG. 1: the marginalized joint likelihood for the lensing amplitude Al and the scalar spectral index Ua panel (a), the amplitude 
of the (linear) power spectrum on the scale of 8/i“^Mpc, erg panel (b) and the reionization optical depth r panel (e). In all 
three panels different colors correspond to different combination of mock data sets as shown in legend. 

The darker and lighter shades correspond respectively to the 68% C.L. and the 95% C.L.. 


TABLE II: Fiducial parameters of the mock data sets 


CP 

Mock-A 

Mock-B 

10'’A, 

2.10745 

2.14338 

Ua 

0.97468 

0.97156 

T 

0.0611 

0.0664 

Q,bh? 

0.022674 

0.022379 

Uch'^ 

0.11639 

0.11748 

Ho 

69.02 

68.39 


1.28 

1.02 


0.06 

0.06 


“Fiducial parameter values in Mock A data sets are the best-fit 
values of base-ACDM-|-AL to Planck-2015 TT-|-lowP data sets. 

^’Fiducial parameter values in Mock B data sets are the best-fit 
values of base-ACDM-t-An to PZancfc-2015 TT-|-lowP-|-lensing data 
sets. 


III. MODIFIED GRAVITY IN THE EFFECTIVE 
FIELD THEORY FRAMEWORK 

The clustering of matter curves the trajectories of 
CMB photons and mixes the CMB photon anisotropies 
from different directions. Since the characteristic scale 
of clustering of galaxies is smaller than 10Mpc/i“^ we 
can use the flat sky approximation to compute the lens¬ 
ing contributions to the CMB anisotropy spectrum. The 
lensing signal, in general, will be encoded in all the types 
of CMB anisotropy modes (T,E,B), and here we only 
show the effect on the TT spectrum as an example, 

cr ^ cr+1 

- cr I ( 4 ) 


in the numerical analysis that follows we take lensing 
effect into account also for TE and EE spectra. For the 
expressions of lensed TE and EE spectra we refer the 
reader to [3]. 

In the one extra parameter extension of the base- 
ACDM model in the Planck-2013 [TU] and Planck-2015 [T] 
results, the collaboration studied the case of varying 
the lensing amplitude parameter, Ap, in the CMB 
anisotropies, which was originally introduced in m- 
This phenomenological parameter is defined by Cf^ —)■ 
Ai^C^, which simply rescales the lensing amplitude con¬ 
tribution to the CMB anisotropies. This parameter, how¬ 
ever, only modulates the CMB anisotropy spectra, Cj^, 
Cj^ and cr! rescales the lensing potential spec¬ 
trum Cf^ but does not rescale the estimator of the lens¬ 
ing spectrum Aff which is computed from the CMB 
anisotropy trispectra [laiTs]. 

In the following we shall elaborate on the relationship 
between this phenomenological parameter and modifica¬ 
tions of gravity. We shall consider two models that are 
enclosed in the background part of the EFT formalism 
for cosmic acceleration [H US- Both can be derived 
from this action written in the unitary gauge and Jordan 
frame 

s = j d'^xyT-g I ^ + rr R + 

-a^c{T)5g°° + ...]+Sm[9r^ (5) 

where we have used conformal time and A and c are 
free functions of time which multiply all the operators 
that are consistent with time-dependent spatial diffeo- 
morphism invariance and contribute to the background 
evolution. The ellipsis indicate operators which would 
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FIG. 2: the effective lensing amplitude A^{1) = C'/‘^(MG)/C/'^(GR) as a function of scale in the modified gravity models 
considered in this paper. Different colors correspond to different values of the model parameters as shown in legend. 


affect only linear and non-linear cosmological perturba¬ 
tions, while Sm indicates the action for all matter fields: 
cold dark matter, baryons, massive and massless neutri¬ 
nos and photons. 

The first modified gravity model that we use is f(R), 
its mapping to the EFT framework was presented in m 
and we refer the reader to [EHll] for detailed discussions 
on the cosmology of these models. 

The second model instead consist in taking a constant 
value for the conformal coupling 17 (a) = and re¬ 

quiring the expansion history to be exactly that of the 
ACDM model. This requirement will then fix, through 
the Friedmann equations, the time dependence of the op¬ 
erators c and A. 

We highlight here that the constant 17 model is not a 
simple redefinition of the gravitational constant. In fact 
the requirement of having a ACDM background with a 
non-vanishing 17, that would change the expansion his¬ 
tory, means that a scalar field is sourced in order to com¬ 
pensate this change. This scalar field will then interact 
with the other matter fields and modify the behaviour of 
cosmological perturbations and consequently the CMB 
power spectra and the growth of structure. For instance, 
it is easy to show that in the constant 17 model, c(r), 
which is vanishing in general relativity, is non-zero and 
reads 

17 

C = -^{Pm + Pm) ■ (6) 

Another general remark we would like to make on the 
models that we consider here, is that they display a rad¬ 
ically different cosmology, as they correspond to two dif¬ 
ferent behaviour of the perturbation’s effective gravita¬ 
tional constant. Viable models, in the f{R) case [l7j, 
correspond to an enhancement of the gravitational con¬ 
stant which in turn results in the amplification of the 


growth of structure that enhances substantially the lens¬ 
ing of the CMB. 

In the second case we consider two possibilities. If the 
constant 17 is positive the model will be characterized by 
a smaller effective gravitational constant resulting in a 
suppression of the growth and consequently a suppres¬ 
sion of the CMB lensing. We shall call this case the 17+ 
model. If the constant 17 is negative, on the other hand, 
the model will have an enhanced effective gravitational 
constant with a phenomenology similar to that of f{R) 
models. In contrast to what happens to the 17+ case, 
that respects all the usual requirements of physical vi¬ 
ability [ini, this model is only classically stable. This 
means that perturbations around the FRW background 
are stable and well behaved but, for example, the sign of 
the scalar field kinetic term is wrong. We shall call the 
case in which the constant 17 can be greater and smaller 
than zero the 17+ model. 

To study the phenomenology of these three models 
we use the EFTCAMB code ^ [201 HI] that allows to 
compute cosmological observables for all the theories en¬ 
closed in the EFT formalism once a precise mapping or 
parametrization is given. An important feature of EFT¬ 
CAMB, besides its versatility, is that it evolves the full 
perturbation equations, on all linear scales, without re¬ 
lying on any quasi-static approximation. In addition it 
checks the stability conditions of perturbations in the 
dark sector in order to ensure that the underlying gravi¬ 
tational theory is acceptable. The implementation details 
of the code and the equations that are solved can be seen 
in |22j with further comments on the models that we are 
considering. 


^ Publicly available at: http://wwwhonie.lorentz.leidenuiiiv.nl/ 
~hu/codes/ 
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a) designer f(R) b) constant Q.± c) constant Q+ 



H Mock-A: TT (Al = 1.28) H Mock-C: XT (Al = 1.28) and lensing (Al = 1.02) 


FIG. 3: the marginalized joint likelihood for the amplitude of the (linear) power spectrum on the scale of 8h“^Mpc, erg and 
the present day value of Log^gBo panel (a), the present day value of the conformal coupling in the case in which it is 

allowed to have positive and negative values panel (b) and in the case in which it is restricted to positive values panel (c). In 
all three panels different colors correspond to different combination of mock data sets as shown in legend and the grey band is 
the marginalized Icr bound on erg from the base-ACDM+An model. The darker and lighter shades correspond respectively to 
the 68% C.L. and the 95% C.L.. 


IV. RESULTS AND DISCUSSIONS 

We firstly check the reliability of our simulations. As 
reported in the Table [n| we use the fiducial Al = 1.28 
for the generation of the Mock-A data set. After going 
through the Markov-Chain Monte-Carlo analysis [53], we 
get the marginalized constraint from the Mock-A data set 
as 

Al = 1.31 ± 0.06 , (68%C.L.) . (7) 

The consistency between the input and output values of 
our simulations can also be seen in the panels of Fig¬ 
ure]^ In addition from panel (a) we can see that there is 
a relatively large positive correlation between Us and Al. 
This happens because a large value of rig will enhance the 
high-t' multipoles, while an increase in Al will smear the 
peaks in the same multipole range. Panel (b) instead 
shows that erg and Al are anti-correlated. An increase of 
Al will suppress the growth inferred from the TT power 
spectrum, hence we end up with lower ug value. Panel 
(c), at last, shows that no significant degeneracy arises 
between Al and the re-ionization optical depth parame¬ 
ter T. 

From all the panels of Figurej^we can see how different 
cosmological parameters reacts to the use of a tension 
data set. The shift of blue and red contour from Mock-A 
and Mock-C data sets shows the tension between them. 
In particular in panel (a) we can see that the Mock-C 
data set gives a value of Al that is in the middle of 
the Mock-A and Mock-B data sets that were used to 
build it. From panel (b) and (c) we can instead see that 
the marginal probability distribution of Ug and r does 
not change signihcantly as the tension in the data set 


TABLE III: The best-fit 


2 

Xbest —fit 

Mock-A 

Mock-C 

base-ACDM 

1967.373 

3314.547 

base-ACDM+AL 

1951.122 

3303.599 

f{R) 

1952.325 

3305.251 


1967.604 

3314.676 


1951.109 

3304.507 


is introduced. This happens because the constraint on 
these two quantities are dominated by the TT and EE 
power spectra. 

After checking the consistency of our simulations we 
move to their interpretation within the modified gravity 
models here considered. 

First of all we check whether these models can re¬ 
produce the input amplitude of lensing. To do so 
we define an effective lensing amplitude as A£®(£) = 
Cf'^{MG)/Cf‘^(GR) and we plot it for several choices 
of the parameters defining our modified gravity models 
in Figure]^ From both panels we can see that the effec¬ 
tive lensing amplitude is generally a function of scale so 
that the accuracy of using a scale independent approxi¬ 
mation for it is limited by the accuracy of observations. 
That is, if observations are not precise enough then the 
scale independent approximation can work while if the 
observations have enough sensitivity we would bias pa¬ 
rameter estimation by considering it scale independent. 
This point will be more and more relevant as the exper¬ 
imental accuracy improves. 

The f{R) case, shown in panel (a), in particular, dis¬ 
plays an Af5 that is monotonically increasing as a func- 
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tion of scale. At low multipoles the agreement with a 
scale independent model is as good as few percent while 
in the high-£ range it is off by 40% depending on the 
present value of the Compton wavelength of the scalaron. 
The constant n_|_ and il,± models are shown in panel (b). 
Unlike the f{R) case, has a wide bump or dip 

around multipoles of few tenths while the scale depen¬ 
dence is somewhat weaker than the previous case. As 
expected, in the U_|_ model the amplitude of lensing is 
decreased as a consequence of the suppression of growth 
while the opposite behaviour shows when Uq is smaller 
than zero. 

We fit these models to our Mock-A and Mock-C data 
sets and we show the marginalized bound on parame¬ 
ters of interest in Figure From panel (a) we can see 
that there is a strong positive correlation between the 
scalaron present Compton wavelength parameter Bq and 
as- This effect is well known and corresponds to the 
fact that as the effective Newton constant is increased 
the growth of matter perturbations is increased as well. 
When considering just Mock-A our results show that in 
order to mimic Al 1.3, the marginalized constraint on 
Bq have to point significantly toward large values of the 
scalaron Compton wavelength. The marginal bound is 
then —1.15 < Logio.So < —0.04 at 95% C.L.. 

Due to the significant degeneracy between Bq and Og 
such values of Bq will lead to strong enhancement of the 
growth rate at relatively small scales making the og value 
too large so that the tension between Planck and and LSS 
surveys, such as CFHTLenS would become even 

worse. For this reason when lensing is added in Mock- 
C, the tension in Al, that for the lensing data set is 
smaller, pushes the posterior of Bq toward smaller values, 
making it closer to GR. A similar effect was also observed 
for the Planck-2015 data set [55] and in the Planck-2013 
data UZHSOj. 

As expected a similar result is found also in the 
model, as shown in panel (b). Here the correlation be¬ 
tween is negative due to the fact that negative val¬ 

ues of U correspond to significant deviations from GR, 
and consequently to stronger enhancement of CMB lens¬ 
ing. All the conclusions, previously discussed in the f{R) 
case, apply also to this model with the relevant exception 
that the enhancement of the growth is not so dramatic 
as in /(i?). The bound on og from the Al fit, shown in 
Figure]^ as a gray band, is in fact almost compatible with 
the one in the case. This is due to the fact that devi¬ 
ations from scale independence of AJ®, in this model, are 
weaker and a constant Q is more efficient at mimicking a 
scale independent Al than f{R) models. Since a weaker 
effective Newton constant is disfavoured by the Planck- 
2015 CMB anisotropy data, the best-fit parameter in the 
model mimic those in the base-ACDM model and the 
model result to be compatible with GR. This conclusion 
is similar to the one found in EIIIMIIIII- 


At last in Table III we list the best-fit for the one 


parameter extensions of the base-ACDM model that we 
investigated in this paper. From there we can see that 
the best-fit values in the base-ACDM-t-Au, f{R), 
models are similar and are noticeably smaller than those 
from base-ACDM or D+ models with ~ —16 from 
the Mock-A and Ay^ ~ —11 from the Mock-C data sets. 
The best-fit y^ in the case is slightly better than 
the f{R) case and closer to the Al one because of the 
weaker scale dependence of A^®. 


In conclusion, Planck-2015 results revealed some inter¬ 
esting tensions between CMB temperature and polariza¬ 
tion anisotropies and CMB lensing. These tensions add 
to the one reported by LSS surveys, like CFHTLenS, that 
seem to favour a smaller og. If this is not due to some 
unaccounted systematic effects, this might be an indica¬ 
tion of exotic physics beyond the base-ACDM model and 
as such should be investigated in great detail. 

In this paper, using a simulation of the Planck-2015 data 
set, we tried to reconcile this tension with some mod¬ 
ified gravity models. In particular we studied whether 
this tension can be mitigated by /(i?) models or models 
characterized by a constant conformal coupling between 
gravitational and matter perturbations. We found that 
the f{R) and the models can mimic the role of Al 
even if, generally, the amplitude of lensing, with respect 
to the GR case, is scale dependent. In order to provide a 
good fit to the Planck-llke CMB anisotropy spectra, how¬ 
ever, these models predict an enhancement of the growth 
on smaller scales that make the tension in tjg even worse. 
From these results we can conclude that the tension be¬ 
tween the growth of matter perturbations assessed from 
CMB power spectra, CMB lensing and LSS surveys can 
be mimicked by modified gravity models but is hardly 
relieved by simple models beyond the standard ACDM 
one. This in turn suggests to investigate more compli¬ 
cated models, possibly with different time dependencies, 
to allow for different regimes of growth at the times at 
which each of these observations is more sensitive. 


Acknowledgments 


We thank Carlo Baccigalupi, Julien Lesgourgues, Mat- 
teo Martinelli, Sabino Matarrese, Laurence Perotto, 
Alessandra Silvestri, Wessel Valkenburg and Matteo Viel 
for helpful correspondences and discussions. B.H. is sup¬ 
ported by the Dutch Foundation for Fundamental Re¬ 
search on Matter (FOM). M.R. acknowledges partial sup¬ 
port from the INFN-INDARK initiative. M.R. thanks 
the Instituut Lorentz (Leiden University) for hospitality 
while this work was being completed. 



7 


[1] P. A. R. Ade et al. [Planck Collaboration], 
arXiv:1502.01589 [astro-ph.COj. 

[2] L. Perotto, J. Lesgourgues, S. Hannestad, H. Tu and 
Y. Y. Y. Wong, JCAP 0610, 013 (2006) [astro- 
ph/0606227]. 

[3] J. Tauber et al. [Planck Collaboration], astro- 
ph/0604069. 

[4] W. Hu and T. Okamoto, Astrophys. J. 574, 566 (2002) 
[astro-ph/0111606]. 

[5] A. Lewis, A. Challinor and A. Lasenby, Astrophys. J. 
538, 473 (2000), ]astro-ph/9911177]. 

[6] http://canib.info . 

[7] M. Tegmark, A. Taylor and A. Heavens, Astrophys. J. 
480, 22 (1997) [astro-ph/9603021]. 

[8] A. Lewis, Phys. Rev. D 71, 083008 (2005) [astro- 
ph/0502469]. 

[9] A. Lewis and A. Challinor, Phys. Rept. 429, 1 (2006) 
[astro-ph/0601594]. 

[10] P. A. R. Ade et al. [Planck Collaboration], Astron. Astro¬ 
phys. 571, A16 (2014) [arXiv: 1303.5076 [astro-ph.CO]]. 

[11] E. Calabrese, A. Slosar, A. Melchiorri, G. F. Smoot 
and O. Zahn, Phys. Rev. D 77, 123531 (2008) 
[arXiv:0803.2309 ]astro-ph]]. 

[12] P. A. R. Ade et al. [Planck Collaboration], 
arXiv:1502.01591 [astro-ph.CO]. 

[13] P. A. R. Ade et al. [Planck Collaboration], Astron. Astro¬ 
phys. 571, A17 (2014) [arXiv: 1303.5077 [astro-ph.CO]]. 

[14] G. Gubitosi, F. Piazza and F. Vernizzi, JGAP 1302, 032 
(2013) ]JCAP 1302, 032 (2013)] [arXiv:1210.0201 [hep- 
th]]. 

[15] J. K. Bloomfield, E. E. Flanagan, M. Park and S. Wat¬ 
son, JCAP 1308, 010 (2013), [arXiv:1211.7054 [astro- 
ph.CO]]. 

[16] Y. S. Song, W. Hu and 1. Sawicki, Phys. Rev. D 75, 
044004 (2007) [astro-ph/0610532]. 

[17] L. Pogosian and A. Silvestri, Phys. Rev. D 77, 
023503 (2008) [Erratum-ibid. D 81, 049901 (2010)] 


[arXiv:0709.0296 ]astro-ph]]. 

[18] R. Bean, D. Bernat, L. Pogosian, A. Silvestri and 
M. Trodden, Phys. Rev. D 75, 064020 (2007) [astro- 
ph/0611321]. 

[19] A. De Felice and S. Tsujikawa, Living Rev. Rel. 13, 3 
(2010) ]arXiv:1002.4928 [gr-qc]]. 

[20] B. Hu, M. Raveri, N. Frusciante and A. Silvestri, Phys. 
Rev. D 89, no. 10, 103530 (2014) [arXiv:1312.5742 [astro- 
ph.CO]]. 

[21] M. Raveri, B. Hu, N. Frusciante and A. Silvestri, Phys. 
Rev. D 90, no. 4, 043513 (2014) [arXiv: 1405.1022 [astro- 
ph.CO]]. 

[22] B. Hu, M. Raveri, N. Frusciante and A. Silvestri, 
arXiv: 1405.3590 [astro-ph.IM]. 

[23] A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002) 
[astro-ph/0205436]. 

[24] C. Heymans, L. Van Waerbeke, L. Miller, T. Erben, 
H. Hildebrandt, H. Hoekstra, T. D. Kitching and Y. Mel- 
lier et al., Mon. Not. Roy. Astron. Soc. 427, 146 (2012) 
[arXiv: 1210.0032 [astro-ph.CO]]. 

[25] T. Erben, H. Hildebrandt, L. Miller, L. van Waerbeke, 
C. Heymans, H. Hoekstra, T. D. Kitching and Y. Mellier 
et al., Mon. Not. Roy. Astron. Soc. 433, 2545 (2013) 
[arXiv:1210.8156 [astro-ph.CO]]. 

[26] P. A. R. Ade et al. [XXX Collaboration], 

arXiv: 1502.01590 [astro-ph.CO]. 

[27] B. Hu, M. Raveri, A. Silvestri and N. Frusciante, 
arXiv: 1410.5807 [astro-ph.CO]. 

[28] J. Dossett, B. Hu and D. Parkinson, JCAP 1403, 046 
(2014) ]arXiv:1401.3980 [astro-ph.CO]]. 

[29] B. Hu, M. Liguori, N. Bartolo and S. Matarrese, Phys. 
Rev. D 88, no. 12, 123514 (2013) [arXiv: 1307.5276 [astro- 
ph.CO]]. 

[30] A. Marchini and V. Salvatelli, Phys. Rev. D 88, no. 2, 
027502 (2013) ]arXiv:1307.2002 [astro-ph.CO]]. 



